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PLATE MOTION AND THE SECULAR SHIFT 
OF THE MEAN POLE 

Han-Shou Liu 
L. Carpenter 
R. W. Agreen 
Geodynamies Branch 

ABSTRACT 

Recent development of the global plate motion indicates 
that changes in the products of inertia of the earth due to 
tectonic plate movement may provide a secular shift of the 
mean pole. In this paper we present a mathematical pro- 
cedure for calculating this shift based on the plate theory. 
Explicit expressions are obtained for the dependence of the 
secular polar shift on the dimensions and locations of the 
plate boundaries. Numerical results show that the secular 
motion of the mean pole is 0"0002 year '' in the direction 
of 67°W. Hence, it is deduced that the influence of the 
plate motion on the secular polar shift may account for 10% 
of the observed value. 
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PLATE MOTION AND THE SECULAR SHIFT OF THE MEAN POLE 


INTRODUCTION 

The analysis of the concurrent observations of the International Latitude 
Service shows that the mean pole of the earth has a secular motion which con- 
sists of a progressive component of about 0"0035 yr _1 along the meridan 65°W. 
(Markowitz, 1968). Subtracting the apparent motion of the mean pole due to the 
local drift, the secular motion of the mean pole is Oi'0022 yr' 1 in the direction 
77°7W (Yumi and Wako, 1968). This mean pole, obtained from a suitable com- 
bination of the polar coordinates covering several years with the periodic com- 
ponents of the annual and Chandlerian motion removed, is the pole of inertia of 
the earth (Melchior, 1972). Recently, Kane (1972) and Bostrom (1971) have 
proposed a link between polar motion and plate tectonics. They have suggested 
that movement of the lithosphere consisting of irregular plates as described by 
Morgan (1968), LePichon (1968), Isacks, et al. (1968), may lead to changes in 
the location of the principal axes of inertia of the earth and hence to drifting of 
the pole. The purpose of the present paper is to analyze the influence of the 
global tectonic plate motion on the shifting of the pole of inertia or the mean 
pole of the earth. 

ANALYSIS 

A general analysis is required to determine the position of the principal 
axes of the earth when the motion of the lithospheric plates is considered. We 
take a reference system oxyz as principal axes of the earth about which its 


1 



moments of inertia are A, B, and C. Let the lithospheric plates undergo a 
small movement, which turns the principal axes through angles £ , £ and V about 
the reference axes and makes the new principal moments of inertia A', B’ 
andC', about the new principal axes x', y' and z'. Because of the plate motion 
we may define the moments and products of inertia about the reference axes as 
A + a, B+b, C+c and E, F, G. Then it is required to find | , £ and rj in terms 
of these quantities. 

We denote i , m and n as the direction cosines of a line through the center 
of the earth with respect to the reference axes. Let 1 + Af, m + Am and n + An 
be the direction cosines of this line with respect to the new principal axes. Then 
the connection between moments and products of inertia with respect to the two 
sets of axes through the center of the earth is governed by (Wittaker, 1944). 

A'(f +Ai ) 2 + B' (m + Am ) 2 + C’ (n +An ) 2 (1) 

= (A+a)i 2 + (B+b) m 2 + (C+c)n 2 - 2Emn- 2 Fni - 2Gfm 

and 

Ai = 77m - £ n 

Am = |n -tU ( 2 ) 

An = £1 - |m 

By substituting equation (2) in equation (1), we obtain 

A' f 2 + B'm 2 + C'n 2 - 2mn (C'-B’) £ -2nf (A'-C') £ - 2fm (B'-A')tj 

(3) 

= ( A+a ) f 2 + (B+b)m 2 + (C+c)n 2 -2Emn- 2Fnf- 2Gim 
to the first order of small quantities. 
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Equation (3) holds for all values of 2 , m, and n. 

Hence 

A' = A+a for 2 = 1 and m = n = 0 

B* = B+b for m = 1 and 2 = n = 0 (4) 

C' = C+c for n = 1 and 2 = m = 0. 


Thus 


£ = 
£ = 
V = 


E 


E 


C’-B* " (C-B) + (c-b) 

F = -F 

A'-C' (C-A) + (c-a) 

G G 


(5) 


B'-A' (B-A) + (b-a) 

For the case B=A, the axes of reference may always be so chosen that G 

is absolutely zero because all axes in the xy plane were originally principal 

axes. Consequently, equation (5) becomes 

E 


1 = 
f = 


(C-A)+(c-a) 

-F 


( 6 ) 


(C-A) + (c-a) 
V = 0 


Therefore the inclination, i, between the new and old principal axes of z' and 
z is 


J = (E 2 +F 2 ) 2 

(C-A) + (c-a) ( 

and the displacement of the new principal axis z' along the meridian measured 
from the plane of xz is 


Q = 


7T + tan 



(8) 
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7 

Equations (7) and (8) express the magnitude and direction of the deflection 
of the polar axis of the earth in terms of changes in the moments and products 
of inertia which may be caused by the tectonic plate motion. 

CHANGES IN INERTIA 

The effect of plate motion in deflecting the pole is equivalent to the move- 
ment of a given quantity of matter from one part on the surface of the earth to 
another. According to the plate theory, the past and present positions of the 
oceanic ridges and trenches and their heights and depths give some idea of the 
amount of matter at disposal. It is required, therefore, to determine the 
magnitude and direction of the possible deflection of the pole which can be caused 
by the movement of given quantities of matter from one part of the earth to an- 
other, subject to the conditions of the height, depth and position of plate 
boundaries. For the present purpose, it is immaterial what may be the forces 
which produce the plate motion; but it is assumed that the plate motion proceeds 
with uniform velocity. 

In the process of plate motion the earth may be considered as composed 
of the original ellipsoid, together with a superposed layer of material which is 
positive in oceanic ridges and negative in oceanic trenches. If we take axes with 
origin at the center of the ellipsoid and symmetrical thereto, and let h$(£,A) 
represent the depth of this layer at point f, A, where £ is the co-latitude and 
A is longitude, the integral of <f>(£,A) over the surface of the ellipsoid should be 
zero because the volume of the earth remains constant. Then by describing 


4 



h<£(£,A) along plate boundaries, the mass redistribution of tectonic plates may 
be expressed. In order to calculate the products and moments of inertia due to 
plate motion, we use spherical coordinates: x = r Sinf CosA, y = r Sinf SinA, 
z = r Cosf, where r is the radius of the earth. An element of volume in these 
coordinates is dV = r 2 Sin£d JdAdr. Therefore the changes in products and 
moments of inertia due to plate motion are 


ft 4 ft 2 

E = hpr 4 / / q>(£,A) Sin 2 £Cos£ SinAd £dA 


-.3 * 1 


ft 4 ft 2 

F = hpr 4 J J $(£,A) Sin 2 £Cos £CosAd£ dA 

" 3 J 1 1 


( 9 ) 


fi 4 fl 2 

a = hpr 4 / J <j>(£,A) Sin£ (Sin 2 £ Sin 2 A +Cos 2 £)d£dA 

t 3 " 1 

fl 4 fi 2 

c = hpr 4 J^ J ^ <J>(£,A) Sin 3 f d£dA 

where p is the density of the surface of the earth, $ (£, A) is a function to deter- 
mine if mass is added by the plate motion or removed and £,, l 2 , £ 3 and t 4 are 
limits on the plate boundaries. 

For the estimation of the secular shift of the mean pole, a computer pro- 
gram is developed to follow the evolution of the products and moments of inertia 
of the earth. $(£,A) may be chosen such that it is zero at the limits of integra- 
tion. In this case a double application of the trapezoidal rule (Hildebrand, 1956) 


gives 
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E = hpr"* 


(N+1)(M+1) 


M N 

*ll[< 

1 = 1 k=l 
M N 


<f>(£ k ,Aj) Sin 2 £ k Cosf k SinXj 


F = hpr ' X 2 [*UitX|) s in i £kCosf k CosX i 


i=i k=i 
M N 


( 10 ) 


IY1 »V r~ 

a = hpr4 ^+lHM^ 2 2 L* (f i«’ x i ) Sin^( Sin2 ^k Sin'Xj +Cos 2 £ k 


i=l k=l 

M N 




i = l k=l 


where 


L = jL &_ ) + 1 


N+l 


1 J 


A +i 
1 M+l 


which are coded for computation by a computer. The formalism outlined here 
can be also applied for the estimation of polar motion due to the inertia changes 
arising over shorter time scales from the accumulation of ground water by 
diverting the flow of large rivers (Batrakov, 1972). 

SHIFT OF THE INERTIA CENTER 

In the derivation of equations (7), (8) and (10), we have assumed that the 
earth’s center of inertia is fixed at the origin. However, during the process of 

plate motion, the coordinates (x,, y,, z ) of the inertia center are shifted by 
X ’ = / / $ (£>X) Sin 2 ? CosA d £ dA 

hpr 3 f C 

Yi = M J J $ (£,A)Sin 2 f SinAdSdA (11) 

z , = ~ “ // $ (t,X)SinfCosXdtdX 
4 

where M = — n p m r 3 is the mass of the earth and p m is the mean density. 
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The products and moments of inertia of the earth about axes parallel to 
x, y, and z through (x Jf y , z x ) in the process of plate motion are 
A' = A+a - M(y, 2 + z, 2 ) = A+a (1- f) 

B' = B4b - M(x, 2 + Zj 2 ) = B+b (1- |) 

C* = C+c - M(x, 2 + y, 2 ) = C+c (l-£) (12) 

E» = E - My, z, = E(l- ~) 

F' - F - Mx,z, = F(l--|) 

Qt Q *y £ £ 

From equations (9), (11) and (12) we note that the ratio of — , — and — 

varies as h/r. The crest of oceanic ridges stands several km higher than the 
flanking abyssal plains (Van Bemmelen, 1972); hence h/r is in the order of 10 ' 3 . 
Therefore, the changes in products and moments of inertia due to the shift of 
the inertia center by plate motion is negligible, and the application of equations 
(7), (8), (9) and (10) for the estimation of the polar motion is justified. 

RESULTS OF COMPUTATIONS 

To make the problem entirely determined we follow the ideas of LePichon 
(1968), Morgan (1968) and Isacks, et al. (1968) to divide the earth's surface into 
plates which stay undeformed except at their boundaries where material may be 
added or moved. This simplification will lead to the mathematical solution 
which can be considered as an approximate solution to the actual problem of the 
redistribution of mass on the surface of the earth due to tectonic plate motion. 
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During the sea-floor spreading process, which occurs at oceanic ridges, 
new crust is generated along plate boundaries. At the deep oceanic trenches, 
plate boundaries are destructive. Chase (1972) and Deffeyes (1970) have calcu- 
lated the rate of crustal flux on the plate boundaries. According to their result, 
a surface area of 5.867 km 2 is involved in the creation and destruction of 
lithosphere each year in connection with plate motion. 

The world-encircling mid-ocean ridge system is a topographic rise, 
thousands of kilometers in width and tens of thousands of kilometers in length. 

If we assume that the net thickness of the lithospheric material which is created 
along the oceanic ridges is 3 km, then the total volume of lithosphere which is 
displaced is 17. 6 km 3 annually. These data were used to compute the cor- 
responding products and moments of inertia as described previously. 

In the process of computation, the following parameters and coefficients 
are used (Allen, 1963) 

A = 0.3295 Mr, 2 , C = 0.3306 Mr, 2 , 

M = 5. 977 x 10 27 g. , r, = 6378.17 km. , 

P = 2. 64 g. cm' 3 , r = 6371. 03 km. 

Taking these data, we have calculated the magnitude and direction of the secular 
shift of the mean pole of the earth by means of equations (7), (8) and (10). The 
results are i = 0!’0002 yr’ 1 and 9 = 67° W. 
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DISCUSSION 


Markowitz (1968) has studied and discussed the reality of the observed 
secular change of the position of the mean pole by assuming that the positions 
of the International Latitude Service stations did not change with each 
other. However, the real change in mean latitude of any station which is ex- 
pected from a supposed movement related to plate motion would produce an 
apparent change of the mean pole (Kaula, 1970). Arur and Mueller (1971) and 
Mueller and Schwarz (1972) have attempted to isolate the observed secular 
motion of the mean pole from plate motion and investigated to what extent the 
secular change is compatible with plate theory. No conclusion can be drawn 
from their study because of the limited available data of large uncertainties. 

The effect of earthquakes on the polar motion (Mansinha and Smylie, 1967) may 
at most account for 30 percent of the observed secular polar shift (Ben-Menahem 
and Isreal, 1970). Moreover, no attempt seems to have been made to explain 
the influence of the site and azimuth of an earthquake on the direction of the 
secular polar, motion; and this point should indeed be considered as an essential 
one (Melchior, 1972). In the present paper we have formulated the effect of 
plate movement on both the magnitude and direction of secular polar motion. 
The changes in the products and moments of inertia of the earth due to plate 
motion are computed for mass redistribution along plate boundaries that have 
been analyzed in the literature. From the changes in the products of inertia, 
the contribution to the resulting polar shift is calculated. Applying the plate 
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theory, we arrived at an annual secular polar shift of 0V0002 in the direction 
67°W. The theoretical foundation of this study is based upon an extremely simple 
rigid global plate model of the earth. In order to make progress in this direc- 
tion, it is indispensable to describe the global tectonic plates in terms of 
rheology (Van Bremmelen, 1972) which could change considerably the order of 
magnitude of the estimated effects. Finally, since the degree of correlation 
between our result and observed secular motion of the mean pole seems to be 
highly dependent on data used, better fundamental observation of polar motion 
with new methods as developed by Smith, et al. (1972) is needed. 
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